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We calculate the energies of three-quark states with definite permutation symmetry (i.e. of SU(6) 
multiplets) in the N=0,l,2 shells, confined by the Y-string three-quark potential. The exact Y- 
string potential consists of one, so-called three-string term, and three angle-dependent two-string 
terms. Due to this technical complication we treat the problem at three increasingly accurate 
levels of approximation: 1) the (approximate) three-string potential expanded to first order in 
trigonometric functions of hyper-spherical angles; 2) the (approximate) three-string potential to all 
orders in the power expansion in hyper-spherical harmonics, but without taking into account the 
transition(s) to two-string potentials; 3) the exact minimal-length string potential to all orders in 
power expansion in hyper-spherical harmonics, and taking into account the transition(s) to two- 
string potentials. We show the general trend of improvement of these approximations: The exact 
non-perturbative corrections to the total energy are of the order of one per cent, as compared with 
approximation 2), yet the exact energy differences between the [20, 1 + ], [70, 2+], [56, 2+], [70, 0+]- 
plets are shifted to 2:2:0.9, from the Bowler and Tynemouth separation rule 2:2:1, which is obeyed 
by approximation 2) at the one per cent level. The precise value of the energy separation of the 
first radial excitation ("Roper") [56', + ]-plet from the [70, l~]-plet depends on the approximation, 
but does not become negative, i.e. the "Roper" remains heavier than the odd-parity [70, l~]-plet in 
all of our approximations. 

PACS numbers: 12.39.Pn,14.20.-c 



I. INTRODUCTION 



The so-called Y-junction string three-quark potential, denned by 
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VV=ermin > Wl~ x o|- (1) 
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i=l 

has long been advertised P, [H, [H, 0, [f| [|| as the natural approximation to the flux tube confinement mechanism, that 
is allegedly active in QCD. Lattice investigations, Refs. UH, @, however, contradict each other in their attempts to 
distinguish between the Y-string, Fig. [TJ and the A-string potential, see Fig. 
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V A = aJ2 l x i- x il: (2) 
i<j=l 

which, in turn, is indistinguishable from the sum of three linear two-body potentials. The present point of view held 
among the lattice QCD practitioners is that there should be a smooth cross-over from the A to the Y-potential at 
interquark distances of around 0.8 fm [irj| |. Exactly how this cross-over should be implemented is not clear just now. 
Moreover, quantum fluctuations of the three flux tubes lead to Luscher type corrections in the potential [TlJ. The 
evaluation of such corrections is beyond the scope of the present paper. 

Over the past 25 years, the Y-string potential has been used in several, more or less systematic studies of baryons 
in the (constituent) quark model with various hyperfine interactions [!, 0, M, HH, and yet some of the most basic 
predictions of this potential acting alone, i.e. without hyperfine interaction, such as the spectrum of the low-lying 
three-quark states remain unknown. The older calculations 0, 0] treat the Y-string potential only in first-order 
perturbation theory, whereas the more recent ones [l2| invoke equivalence with A string potential 26], up to an 
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FIG. 1: Three-quark Y-junction string potential. 




FIG. 2: Three-quark A-shape string potential. 



overall multiplicative factor / ~ 0.5493 (| < / < ^=) for the string tension a. That leads us to a third set of 
calculations that rely on the (as yet not proven) "equivalence" of the Y-string_potential and the A string one to look 
for the "best two-body approximants" to the Y-string potential, see Refs. |13l . |22T ] and follow-up references. Note, 
however, that these papers use only the hyper-radial approximation and they are only concerned with the overall 
strength of the coupling. 

We wish to resolve this dilemma using analytical methods, whenever possible [271 ]. So, at least for the time being, 
we have to resort to calculations in the nonrelativistic quark model. We shall not repeat the two-body A string 
potential calculation, as there is a large body of literature on the subject: instead, we shall use the paper by J.M. 
Richard and P. Taxil fl9j | as a benchmark. That paper uses the hyper-spherical harmonics method, which we also 
use, with only a slightly different practical implementation: whereas Ref. [lj| uses the Jacobi coordinates (see Sect. 
HI)) in the evaluation of the relevant matrix elements, we shall evaluate them directly in the hyperspherical variables. 

For the sake of completeness and clarity we shall start with two extensive technical preparatory sections: First, in 
Sect. |n]we outline the three main technical problems, and then in Sect. IIHI we address them one after another; for 
this reason we divide Sect. Hill in which we formulate the methods that we use, into three parts: the first one, Sect. 
IIII Al is about the angular momentum recoupling algebra necessary to deal with the non-conserved "partial" angular 
momenta, the second part, Sect. IIII B[ deals with the square root(s) in the Y-string potential, and finally the third 
one, Sect. IIII C[ addresses all four forms of the string potential together. Our results are shown in Sect. IIVI which 
is divided into four parts: one part that comprises Sect. IIV Al deals with the first-order perturbative approximation 
results, wherein the angular dependence is treated as a (small) perturbation to the linear potential oscillator, Sect. 
IIV A H and two with non-perturbative ones in Sect. IIVBI one for the Y-string potential in Sect. IIVB II and another 
for the complete string potential in Sect. IIV B 21 The final Section [V] contains a summary and discussion of our 
results. A number of technical issues are discussed in Appendices lAIBICl 
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II. THE Y-STRING POTENTIAL 



The complexity of the potential Eq. ([T]) is best seen when expressed in terms of three-body Jacobi (relative) 
coordinates p, A 



P = ^|( X l ~ X 2)' 

A = -^(xi +x 2 - 2x 3 ), 
as follows. The exact string potential Eq. |T]) consists of the so-called Y-string, or three-string term, 



String = V Y = a x j-(p 2 + X 2 + 2\px A|), 



(3) 
(4) 

(5) 



2p 2 - V3p ■ A > -py V + 3A 2 - 2V3p ■ A 
when { 2p 2 + V3p ■ A > -p\/p 2 + 3A 2 + 2y/3p ■ A 
3A 2 - p 2 > -i v /(p 2 +3A 2 ) 2 -12(p-A) 2 . 

and three angle-dependent two-part string, or the so-called V-string, terms, 



V, 



string 



a [\j\{p 2 + 3A 2 + 2\/3p -\) + \j\{p 2 + 3A 2 - 2V3p ■ A) 
2p 2 - V3p ■ A > -py/p 2 + 3A 2 - 2\/3> ■ A 



(6a) 



string 



string 



when < 



<r\V2p- 



when < 



when < 



2p 2 + V3p ■ A > -py/p 2 + 3A 2 + 2V3p ■ A 
3A 2 - p 2 < —\\Jip" 1 + 3A 2 ) 2 - 12(p • A) 2 

yi(p 2 + 3A 2 + 2V3p.A)^) 

2p 2 - VSp ■ A > -pV p 2 + 3A 2 - 2V3p ■ A 
2p 2 + V3p ■ A < -pVp 2 + 3A 2 + 2V3p ■ A 
3A 2 - p 2 > —\\Jip" 1 + 3A 2 ) 2 — 12(p • A) 2 

yi(p 2 + 3A 2 -2V3p.A)^) 

2p 2 - V3p • A < -p^ /p 2 + 3A 2 - 2V3p ■ A 
2p 2 + V3p ■ A > -py/p 2 + 3A 2 + 2V3p ■ A 
3A 2 - p 2 > -iv/(p 2 + 3A 2 ) 2 -12(p-A) 2 



(6b) 



(6c) 



Here, the reasons for the lack of use of the exact potential Eq. ([T]) become clear: i) it is a genuine three-body operator 
with a complicated and unusual ( "area term" ) angular dependence under the square-root of the most important term 
(the Y-junction string potential) that leads to the non-conservation of the individual Jacobi coordinates' angular 
momenta [14( and hugely complicates the equations of motion; ii) the square-roots appearing in all four functional 
forms of the potential make this task even more difficult; iii) the presence of four different functional forms of the 
potential depending on the configuration space angles makes the integration of the equations of motion difficult as 
one cannot easily separate the angular and radial integrals. 

Perhaps the simplest, yet realistic approximation to the exact string potential Eq. |T]) is the Y-string, or the 
"three-string" potential, Eq. ([5]), that is used in the whole configuration space, i.e., even when one of the angles in 
the triangle exceeds 120°. In that way one avoids the cumbersome transition to the V-string potentials, see problem 
iii) above [28]. Still, even this simplified approximation suffers from two difficulties mentioned above: i) an unusual 
( "area term" ) angular dependence under the square-root that leads to the non-conservation of the individual Jacobi 
coordinates' angular momenta; ii) the square-root. We shall address these problems in successive steps: i) the area 
term turns out to be exactly (analytically) integrable, but it requires complicated angular momenta recoupling algebra 
and the exact value of a particular one dimensional angular integral. Problem ii), the square root, can be treated, at 
first, by a series expansion, i.e. in perturbation theory, and then by numerical evaluation of the complete functional 
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expression, i.e. in non-perturbative approximation. Finally, the last issue iii) is tackled, at the price of considerable 
inconvenience: the necessary hyperangular matrix elements of the complete three-body potential can be evaluated 
using above mentioned methods, except that two previously separate integrals (over one quasi-radial and one angular 
variable) have to be performed simultaneously, as the integration boundaries involve both variables. 

It turns out that the crucial ingredient for the success of this effort is the application of the so-called hyper-spherical 
coordinates/angles |l7} . or a particular variation thereof, the cosines of the relative angle 9 between the Jacobi 
coordinates p, A and of the angle 1\ defined by way of the ratio of the moduli p, A of the two Jacobi coordinates. 

The hyper-spherical method has been used widely in the few-body atomic and nuclear physics, for review see Ref. 
[24j . but we are aware of only one paper, Ref. [22j], that uses it in the context of the three-quark Y-string problem, 
which is perhaps ironic, as it turns out to be the most natural set of coordinates for the problem at hand, whereas 
its use in two-body atomic and nuclear potential problems is plagued by slow con verg ence of the expansion. Some 
papers have applied the hyper-radial approximation to the three-quark problem [l3l Il9j | but that misses the essential 
points discussed here. 



III. METHODS 



The first problem (angular momentum) is generic to all (string) three-body potentials and is thus independent 
of the approximation used, so its solution will be used subsequently in both perturbative and non-perturbative 
approximations. Moreover, its solution can be incorporated into the hyper-spherical formalism that will used later. 
For this reason we start with the angular momentum recoupling. 



A. Angular momentum recoupling 

The 'vector cross-product, or the "area term" 2 \p x A | in this potential has some curious properties in classical 
and quantum mechanics: it conserves the sum L = l p + l p of the two partial (orbital) angular momenta l p — p x p p 
and 1a = A x pa, but not their difference flij . i.e., the individual (orbital) angular momenta are not conserved. As a 
consequence, there can be "spilling" of the orbital angular momentum from the p (normal) mode into the A one and 
vice versa. Only the radial excitations of the S-wave ground state remain immune to this spillage. 

We need to evaluate matrix elements of the following form 

l 2f ) L f , M f \\ sine\\ (l u ®l 2i )L h Mi) = ^ {l lf ,m lf ;l 2f ,m 2f \L f ,M f ) 

(lii,mii;l 2 i,m 2 i\Li,Mi) x (lif,mif,l 2 f,m 2 f\\ WL6\\lu,m u \l 2 i,m 2 i) (7) 

where 



and 



|pxA| = |pX A| = I sin 01 = \/l - cos 2 i 
pA 



cos0 = cos^i cos6 l 2 + cos((/3i — ip 2 ) sin#i sin9 2 



((h f ® l 2f )Lf,M f \\wa.e\\(l li ® l 2i )Li,Mi) = 5 LfM 8 Mf ,M> ^ v + l ^ 2h f + W 2 *" + l )( 2 ^ + 1) 



L 



(8) 



y j [Y 1 {9 1 ,<p 1 )®Y 1 {9 2l <t> 2 )] L=a (9) 
that is to be inserted into 

(hf,mif,l 2 f, m 2 f\\ sin0||(Zii, run, l 2 i, m 2i )) = 
dn x J dna^^^i^Oyiy^Cfl!,^)^^/^.^^.^^,^) x |sin0| (10) 
Capstick and Isgur, Ref. [J] have reduced the angular matrix elements Eq. (O to 



x (_i)£+*i/+*2, w{ i Uj Lj Ui hf]hf) l 2l )(l lfi o, l lh Q\L, 0)(Z 2/ , 0, l 2u Q\L, 0) 

1 f 1 

x - / d cos 9P L (cos 9) | sin 9\ (11) 

2 J-x 
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a sum of products of SU(2) Clebsch-Gordan (hf, 0, In, 0\L, 0)(hf, 0, hi, 0|L, 0) and Racah Wjlu.i L, Lj hf', h f, hi) 
coefficients, which can be found in standard collections of angular momentum tables, such as Ref. [201 ] . and an 
integral over even-L order Legendre polynomials Pl{x) that can be reduced to a ratio of double- factorial functions 
HI: 



1 A 

- / d cos 8P L (cos 6) VI - cos 2 = 

2 J-i 







for 
for 



odd L : 
L = 



7r (L-l)!!(£-3)!! 
2L (L+2)!!(L-2)!! 



1,3,5,. 



for even L = 2,4,6, ... ^0 



(12) 



In Table U we show the dependence of these matrix elements on the "partial" orbital angular momenta l p ,l\ - note 
in particular the non- vanishing off-diagonal matrix element that ensures the existence of the mass-splitting "mixing" 
terms in the [56, 2 + ] and the [70, 2 + ]-plets. Note that all of the angular matrix elements are exact integer or fractional 
multiples of the basic unit (j), thus leading to the analytic form of the angular matrix elements, shown in the eighth 
column of Table Q] 



TABLE I: Non- vanishing diagonal and off-diagonal matrix elements (| sin#|) = {(hf ® hf) 
P- and D-wave states, and partial waves hf = lif = hi = hi = 0, 1, 2. 



N 



Lf 



IpfX hf Li hi x h 



'A/ 



0x0 



x 1 
x 
x 1 
x 1 
x 
x 



0x0 
1 x 1 
1 x 
1 x 1 
1 x 1 
2x0 
0x2 



1 

3 
4 
1 
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J) 
40 
1 
1 

_a_ 



l\\(lu<S>hi)r), in total S- 



B. Perturbative approximation 

We must first determine the qualitative features of the "area term" 2 \ p x A | and its angular momentum dependence 
on the ordering of states in quantum mechanics, e.g. whether the S-wave states, such as the Roper, are lowered or 
raised in energy as compared with other states in the N=2 band, e.g. the P- and D-waves? 

A historical remark seems in order now: It has been known at least since the late 1970's, see Refs. [HI, [l6l . 
that arbitrary an-harmonic (both two- and three-body) potentials split the N=2 shell harmonic oscillator states |29j 
according to their spatial permutation symmetry classes, or, what is the same, according to the SU(6) multiplets. This 
splitting of the harmonic oscillator spectrum has been worked out to various degrees of mathematical sophistication 
m Refs. [lljllllini for arbitrary two-body potentials, but has been only briefly mentioned in the case of three- 
body ones [l8|. In particular Bowler and Tynemouth [18| have shown, on the basis of the Sp(12,R) group theory 
and first-order perturbation theory, that the energy ordering and splitting of four of the five SU(6) multiplets in the 
N=2 band always remain the same for arbitrary permutation-symmetric three-body potentials, the only exception 
being the energy of the [56, + ] multiplet (containing the Roper resonance), which remains unconstrained by this 
theorem. In other words, the key to the mystery of the Roper's abnormally low mass may well reside in the form of 
the three-quark potential. 



1. The square root in the Y-string potential 



Perhaps the simplest way to address these questions is to expand the root in the Y-string potential Eq. (f5|) in a 
power series 

= ^ + „( 1 + j£^l + ...). (13) 

and then, keep only a few lowest-order non-trivial terms, apply the first-order perturbation theory. The unperturbed 
potential in Eq. (|13[) is slightly more complicated than the harmonic oscillator, so we break this into two steps: 
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We shall use the quasi-linear hyper-radial potential to define the unperturbed Hamiltonian, which leads us to 
introduce the "hyper-spherical (H) string length" Zh, or the "hyper-radius" R, as follows 

R= \J p 2 + \ 2 = Z H - (14) 

The two Jacobi vectors p, A defined in Eqs. ([3]),(j4|), are shown in Fig. 03 Instead of two three- vectors p and A, 
the hyper-spherical formalism introduces the hyper-radius R and the (new) hyper-angle \ by way of the "polar 
transformation" 



p = Rsinx, A = i?cosx with < \ < 7r /2- (15) 



Thus we see that the expansion in Eq. (|13|) corresponds precisely to the expansion in powers of trigonometric 
functions of hyper-angles: 



V Y = (7^2^(1 + sin 2x| sin 6>|), 

= CT y|i?fl+~sin2x|sin0|+---V (16) 

the so-called hyper-spherical expansion, the first term being the "hyper-radial" term and the rest corresponding to 
higher-order "hyper-spherical harmonics" . 

The excited three-quark states then naturally fall into (non-degenerate) multiplets (labelled by the "grand angular 
momentum" quantum number K) of the broken 0(6) symmetry, whose degeneracy is broken by the hyper-angular 
dependent perturbation, see Appendix [AJ These members/subsets of the (broken) 0(6) multiplets coincide with the 
"ordinary" SU(6) multiplets of standard quark model in the two lowest lying bands (usually denoted by the harmonic 
oscillator label N=0,1), the only distinction being that in the N=2 band the first radial excitation of the N=0 ground 
state, a.k.a. the "Roper", does not fall into the K=2 band, but rather belongs to the K=0, hyper-radial excitation 
Wk = 2 "shell" consisting of exactly one state, see Appendix IA 21 Thus, the symmetrized hyper-spherical harmonics 
are the most natural tool to describe the non-relativistic three-quark SUfs(6) multiplets' wave functions and to study 
their mass splittings. 

So, we shall solve the Schrodinger equation with the linear "hyper-radial" potential, see Appendix IC 1[ use its 
solutions as unperturbed states, and treat the second term in Eq. (|16p as the (lowest-order) perturbation in "hyper- 
spherical harmonics", see Sects. liVAlllVl and Appendix IB 11 



C. Non-perturbative approximations 

The hyper-spherical formalism is often advertised as an economical way to tackle the general three-body problems 
in atomic and nuclear physics, but, in fact, it seems as if it had been tailor-made for the Y-string potential. We refer 
to specialized papers for technical aspects of this method (see for instance Refs. [HI). Here we recall just what is 
needed for our purposes; then, we use it below to tackle the problem non-perturbatively. 
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1. V '-string problem in hyper- spherical coordinates 



We expand the three-quark wave function in hyper-spherical coordinates as 

*(R, X ,p,\) = ^VcWc(x,P,A) 



(17) 



where c = Zi, Z2, Jn and if = l\ + l 2 + 2m (occasionally one uses [K] to denote the complete set of hyperspherical 
quantum numbers K, ( 2 ,L, to) and 



y c (x,p = n p ,\ = n x ) = u lul2 , m { x )[Y h {p)®Y l2 (k)] L 

3 

u h ,h,m(x) = N hihim F(-m, h+l 2 +m + 2\l 2 + -| sin 2 x)(cosx) h (sinx)' 2 



I 2(2to + h + l 2 + 2)r(m + h + l 2 + 2)r(rn + l 2 + 3/2) 



r(/ 2 + 3/2) 

where F(—n, n + a\c\z) is the Jacobi function 



F(-n,n + a| C |s) = 



m\T(m + h +3/2) 
n!r(c)r(a + n + s) 



s=0 



(n- s)!s!r(a + n)r(c + s) 



(18) 
(19) 

(20) 



(21) 



and F(c) is Euler's Gamma function. 

An important property of the hyper-spherical formalism is that a complicated three-body problem reduces to a 
set of (coupled) differential equations involving only the hyper-radial matrix elements of the potential V{R)[k][k'] = 
(^(x,0)^ K '^\V(R,x,0)\^(x,0)[K}}hy P ~ang.- The Schrodinger equation of the three-quark system then becomes a set 
of infinitely many coupled equations, 



1 

271 



d 2 1 5 d K(K + 4) _ 
dR 2 RdR R 2 



MR) 



<jRj2Cc,c>ipc>{R) = 



(22) 



with 



Cc,c = 0)c I a/1 + sin 2x1 sin | |tf (x, 0) c ) hyp-ang. 

/ >7r / 2 /■ 

= J cos x 2 sin x 2 c?X j dn p dQ x y*, y/l + sin 2%| sin 0i 2 



(23) 



where cos #12 = p ■ A. 

In special cases, such as the present one, some (finite) subset(s) of equations may decouple, due to the symmetries 
of the interaction potential. The spectrum of the system is then reduced to finding the eigenvalues of the following 
differential equation 



5 d K(K + 4) 



dR? R dR 



- a c R + k 2 



where 



R 2 



a c = 2p\j -crCc.c 



MR) = 



2pE 



(24) 

(25) 
(26) 



and the coupling matrix C c 



r-Tl/2 f-1 

C c ',c = ^Yl Xc '^ L ' J u *c'(x)Mx) cos 2 X sin 2 xd X J dcos6 12 P L <(cos6 12 )yfl + sin2 X \ sin9 12 \ (27) 



can be evaluated using the by now familiar angular re-coupling matrix X c > iCi li 



X d ^ v = y/{2l v + l)(2h + l)(2h + l)(2/ 2 + l)/2 

x (-l) L '+ l i'+ 1 *' W{h, L', L, h;lv,l 2 ){lv,0, h,0\L\ 0)(h,0, l 2 , 0\L', 0) 



(28) 
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FIG. 4: The boundary in the x vs - plane, between the regions in which the two- and the three-string potentials are appropriate, 
see Eqs. J3DJ. 



and the following two-dimensional integrals 

/ P L (cosd)dcos9 f 2 \u c (x)\ 2 \/l + \sm6\sm 2 X (cos* sinx) 2 d X (29) 
J -l Jo 

We may numerically evaluate the necessary integrals, see Appendix IB 21 so as to evaluate the hyper-angular matrix 
elements and recast the Schrodinger equation with the exact potential Vy into a set of (coupled) hyper-radial equations 
for each "hyper-spherical harmonic" that can be solved numerically, see Appendix [Cj 



2. The complete string potential in hyper-spherical coordinates 



To finally solve the exact string problem, one must relax the Y-string potential approximation Eq. ([5]) and use the 
exact potential, i.e., all four of its incarnations, Eqs. (|5l6al6bl6c[) depending on the angles formed by the three quarks. 
This means that one must first determine the boundary in the x vs - plane between the regions in which the two- 
and the three-string potentials are appropriate, see Eqs. (|30[) . There are three such boundaries, determined by the 
three inequalities, that merge continuously one into another at two "contact points" , see Fig. 2] 



cotxi(0) = 
cot% 2 (0) = 



-1 



•\/3 cos 9 + sin ' 
1 



V3 



cos( 



cot X 3(0) = -\Jh~2 cos 2 6 - 2 | sin 9\ V '4 - cos 2 9. (30) 
3 



The hyper-angular matrix elements with these boundary conditions are evaluated in Appendix IB 31 The rest of the 
calculation proceeds as in the case of the Y-string. 



IV. RESULTS 



A. The first-order perturbation theory 



1. Y-string three-quark potential expanded to first powers of hyper-angles 

In Appendix HIT C 21 we have calculated the matrix elements (N K , K [SU FS (6), L p ] \ H££±[ | N K , K [S , C/ FS (6), L p ] > 
for the states/ SU FS (6) multiplets of the three lowest N K (=0,1,2) bands: | N K = 0,K = 0, [56,0+]), |N K = 1,K = 
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1, [70, 1"] ), and |N K = 1,K = 0,[56',0+]), | N K = 2, K = 2, [56,2+]), |N K - 2,K = 2, [70,0+]), |N K = 2,K = 

2, [70, 2+] ), | Nk = 2, K = 2, [20, 1+] ) where we separate each matrix element into its hyper-radial and hyper-angular 
parts: 

(^N K ,[K]I ^ IV'NkJK]) = (#j k ,k|-R|V'N 

k,K )hyp— rad ( ^[K] I sin2x sin6>p A | V>[K] ) hyp-ang, (31) 

where [K] denotes all the quantum numbers, such as the L,l p ,l\ and their magnetic quantum numbers, associated 
with K. The hyper-radial matrix element can be calculated using the virial theorem, see Appendix IC 21 , as follows: 



/3 2 2 

2 ( V'Nk.K I R I V'Nk.K )hyp-rad = g ( V'Nk.K I # I V'Nk.K )hyp-rad = g-^NK.K- (32) 

The calculated energies of states with various values of K and L are listed in Table [TT] Note that the hyper-radial 



TABLE II: The values of the unperturbed energy -En k ,k, the total energy Ek + SEk.l (in units of |-^LJ ) and the 

three-body potential matrix elements' hyper-angular perturbation ( 2 ^* X ^ )h yP -ang. = (sin2x | sin#[)h yp -ang. as well as some 
intermediate steps, for the various K = 0,1,2 states (with all allowed orbital waves L). 



K N K 


E m 


[SU(6),L P ] 


| (sin 2\ | sin ( 


?l> <i,K+^N K ,K,L 


= -En k ,k,l 





3.8175 


[56, + ] 


'2 

li 
■} 

1 


y 3.8175 


4.6658 


1 


4.6582 


[70, 1-] 


S 4.6582 


5.6934 


1 


5.2630 


[56,0+] 


2. 5.2630 


6.4326 


2 


5.4290 


[70,0+] 




|f 5.4290 


6.3942 


2 


5.4290 


[56,2+] 


15 


HI 5.4290 
U 5.4290 
i| 5.4290 


6.4907 


2 
2 


5.4290 
5.4290 


[70,2+] 
[20, 1+] 


W 

225 
1K 


6.6837 
6.8767 



matrix elements of the linear hyper-radial potential are identical for all K=2 band hyper-radial (ground) states. 
Therefore, all the energy differences among various K=2 multiplets are integer multiples of the energy splitting 

"unit" Ay = ^ay / '|(i?)N j K = 0.0327er(-R)N,K = 2§5^n k .k> which is just another manifestation of the validity of the 
Bowler- Tynemouth theorem [18J for this three-body potential solved in hyper-spherical coordinates. 

This theorem has already been confirmed by J.M. Richard and P. Taxil [l9|, in the hyper-spherical formalism with 
linear two-body potentials. Our new contribution here is the (first) proof of this theorem, with the "area term" 
three-body potential as a perturbation, that holds even with non- harmonic oscillator unperturbed states. As in the 
case of Bowler and Tynemouth [l8| the crucial ingredient of the proof is the permutation symmetry of the three-body 
potential. Other simple three body potentials that lack this permutation symmetry do not conform to this theorem. 

Returning to the details of the spectrum, note, however, that the Roper [56, 0+] multiplet (6.433) is roughly half-way 
between the [70, 0+] (6.394) and the [56, 2+] (6.491) multiplets, and for most practical intents and purposes degenerate 
with them both, as their mass difference is just one Ay-. The ratios of the K=2 and other K- values hyper-radial matrix 
elements need not be integers, or rational numbers any more, except in special cases, like in the harmonic oscillator, 
or in the Coulomb potential. In the Y-string potential the upward shift of the Roper [56, 0+] (5.263 — > 6.433) is 
proportional to its unperturbed energy (4.082), as is the upward shift of the K=l odd-parity states (4.658 — > 5.693), 
and the coefficient of proportionality ^ is equal for these two states. 

In other words, the Y-string potential does not move the Roper below the K—l odd-parity resonances, at least 
not in perturbation theory. Note, however, that the largest correction in this hyper-spherical first-order perturbative 
approximation relative to the unperturbed value is 27%, which still does not justify a perturbative treatment. We 
shall therefore try and estimate the effects of the exact potential, i.e. the whole power series at once. 

B. Non-perturbative results 

1. The three-string potential 



Instead of expanding \A + i, where x = sin 6 sin 2\, in a power series we may numerically integrate the double 
integral 

/ P L (cosd) sin Ode [* |u c ,k,l(x)| 2 V 1 + sin d sin 2 X (cosx sinx) 2 d X - (33) 
Jo Jo 
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Such hyper-angular potential matrix elements are the coefficients multiplying the linear hyper-radial potential that 
appears in the (new) hyper-radial Schrodinger equation; that equation, in turn, can be solved exactly, i.e. with- 
out the use of the perturbation theory, and the resulting energy eigenvalues are listed in Table IIIII We can com- 



TABLE III: The values of the unperturbed energy E^ K , the total energy Ek (in units of |-^j=j ) and the three-body 
potential matrix elements' hyper-angular exact non-perturbative matrix element {VV)hyp— ang. 

= (^/l + sin 2 X | sin 0|) 

hyp — ang. j 

for the various K = 0,1,2 states (with all allowed orbital waves L). 



K 


N K 




[SU(6),L P ] 


( sj\ + sin 2x | sin | ) h yP -ang 




= £n k ,k,l 








3.8175 


[56, + ] 


1.2876 


1.18355x3.8175 


4.5182 


1 





4.6582 


[70, 1"] 


1.2876 


1.18355x4.6582 


5.5132 





1 


5.2630 


[56,0+] 


1.2876 


1.18355x5.2630 


6.2290 


2 





5.4290 


[70, 0+] 


1.2350 


1.15110x5.4290 


6.2493 


2 





5.4290 


[56, 2+] 


1.2560 


1.16411x5.4290 


6.3199 


2 





5.4290 


[70, 2+] 


1.2981 


1.18998x5.4290 


6.4604 


2 





5.4290 


[20, 1+] 


1.3402 


1.21557x5.4290 


6.5993 



pare these exact results with the first-order perturbation theory: the energy eigenvalues, after rescaling, see Ap- 

2 

pcndix lUTj are proportional to E ~ (yl ~^7^W^ + sin 2%| sin0|)h yp -angj (modulo an overall hyper-radial de- 
pendent factor that is the same for all K=2 ground states). This can be (double) Taylor-expanded and yields 

2 2 

(V^T/It (! + 5( sin2 xl sin0|) h yp-ang.)) " ^ (y/f 7^=) 3 (1 + i<sin2x| sin6>|) h yp-ang.) , which is the first perturbative 
correction. 

Thus we may expect the exact result to be smaller than the first-order perturbative result on two accounts: 1) from 
the inclusion of all orders in the expansion of the square-root in the potential: the even- and odd-order terms have 
opposite signs (this is an alternating series), see e.g. the first ten terms 

, x x 2 x 3 5 a; 4 7 x 5 21 x 6 33 x 7 429 a; 8 715 a; 9 ml 10 , 

Vl + x = H 1 1 1 1 \-0(x 10 ), 34) 

2 8 16 128 256 1024 2048 32768 65536 y ' v ' 

2 

with rapidly decaying coefficients, meaning that the most important correction to the linear ^ term is the negative ^- 
one; and 2) from the inclusion of all orders in the expansion of the two-thirds-root. Manifestly, both of these effects 
lead to the breaking of the "integer-value rule" of Bowler theorem for the energy splittings of the K=2 multiplets, 
but only at the 1% level. 

Note the overall reduction of the potential matrix elements, as compared with the first-order perturbation theory 
results: the largest correction in the hyper-spherical non-perturbative approach relative to the unperturbed value is 
22%, which is less than the 27% in the hyper-spherical non-perturbative approximation. This fact justifies ex post 
facto the perturbative treatment in Sect. IIV Al 



2. The complete string results 

Finally, the results of the solution to the complete string potential problem are shown in Table IIV1 which shows a 
slight 0(< 1%), yet clear increase, across the board, of the total energy over the Y-string approximation results in 
Table In] 

The overall shifts of energies in this calculation would hardly justify the effort it took to complete it, were it not for 
their (relative) effects on the K=2 level splittings, which are substantial: The violations of the "Bowler- Tynemouth 
theorem" increase to 13% and 9% for the [20, l + ]-[70, 2 + ], and [70, 2 + ]-[56, 2 + ] mass differences, respectively, as 
compared with the [56, 2+]-[70, 0+] one. Note, that the Roper multiplet [56, 0+] (6.2340) has moved below the [70, 0+] 
(6.2665) multiplet, but remains, for all practical purposes degenerate with it. 



V. SUMMARY AND DISCUSSION 



In summary, we have studied the low-lying states in the three quark spectra confined by a pure Y-string potential, 
i.e. without any two-quark potentials, in three different approximations, see Table fVl and Fig. 
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TABLE IV: The values of the unperturbed energy E^ K , the total energy Ek (in units of |-^L j ) and the exact three-body 

potential hyper-angular non- perturb at ive matrix element (Vstring)hyp— ang. ; taken from Table [Xj as well as some intermediate 
steps, for the various K = 0,1,2 states (with all allowed orbital waves L). 



K 


N K 


^Nk.K 


[SU(6),L P ] 


{ V^tring ) hyp — ang 


\ •'string / -C'k 


= £n k ,k,l 








3.8175 


[56, 0+] 


1.2891 


1.18449x3.8175 


4.5218 


1 





4.6582 


[70, 1"] 


1.2891 


1.18449x4.6582 


5.5176 





1 


5.2630 


[56,0+] 


1.2891 


1.18449x5.2630 


6.2340 


2 





5.4290 


[70, 0+] 


1.2401 


1.15943x5.4290 


6.2665 


2 





5.4290 


[56, 2+] 


1.2584 


1.16886x5.4290 


6.3279 


2 





5.4290 


[70, 2+] 


1.2985 


1.19022x5.4290 


6.4617 


2 





5.4290 


[20, 1+] 


1.3404 


1.21567x5.4290 


6.5999 



An attractive Y-string potential splits the N=K=2 band states into degenerate SUfs(6) multiplets: 
[20, 1 + ], [70, 2+], [56, 2+], [70, + ], ordered in descending mass, and following approximately the separation rule of 
2:2:1. The best accuracy results lead to 2.25:2.18:1 splitting i.e. to a maximum violation of this rule is less than 13%. 

The mass difference between the first (hyper-) radial excitation of the ground state, that is the "Roper multiplet" 
[56 ,0 + ], and the odd-parity K=N=1 [70, 1~] multiplet is entirely determined by the difference between and the first 
(hyper-) angular and the first (hyper-) radial excitation eigen-energies in a linearly rising hyper-radial potential, which 
is always negative. In other words, the Roper resonance cannot be lowered below the odd-parity K=N=1 states in 
this potential, irrespective of the string tension constant and the quark masses, which are the only free parameters. 



TABLE V: The eigen-energies (in units of |-^j=j ) of the unperturbed solution to the hyper-central approximation, (see 
the text) K , one perturbative (E^ K L ), and two non-perturbative approximations K L , E^ K L ) to the Y-string 

potential, where the fourth one E^ ' is the exact (numerical) result, for the various low-lying K = 0,1,2 states (with all allowed 
orbital waves L). The last column shows the results for the A-string potential of J.M. Richard and Taxil [l9| . 



K N K 


[SU(6),L P ] 


E {0> 


% Kl K,L 


C N K ,K,L 


E (Y) 










[56,0+] 


3.8175 


4.6658 


4.5182 


4.5218 


5.3592 


1 





[70,1-] 


4.6582 


5.6934 


5.5132 


5.5176 


6.5395 





1 


[56,0+] 


5.2630 


6.4326 


6.2290 


6.2340 


7.3885 


2 





[70,0+] 


5.4290 


6.3942 


6.2493 


6.2665 


7.5409 


2 





[56,2+] 


5.4290 


6.4907 


6.3199 


6.3279 


7.5731 


2 





[70,2+] 


5.4290 


6.6837 


6.4604 


6.4617 


7.6377 


2 





[20, 1+] 


5.4290 


6.8767 


6.5993 


6.5999 


7.7022 



Next we compare our Y-string three-quark potential results with the A-shaped string ones, as obtained in Ref. 
[l9j |. A quick look at Table fVl shows that the A-string eigen-energies are substantially higher than the corresponding 
Y-string ones, with identical string tension a . 

In order to check the (re)scaling hypothesis viz. if the low-lying spectra are identical after a (re)scaling of the string 
tension ay — > /ca, we tabulate the rescaled A-string energies next to the exact Y-string ones in Table fVTl The string 

tension scaling factor / is fixed at ^rj^/|l.2891 = 0.775001, which ensures that the three lowest-lying states/bands, 

viz. (K = 0, Nk = 0, 1) and (K = 1, Nk = 0), have identical eigen-energies, see Table IVTl Note that this value of / 
is substantially different from 0.5493, which was used in Ref. [l2[. The energy rescaling factor is just the | power of 

2 

/: / 2 / 3 = ^7ttj^/|1.2891^ = 0.843726. One can see that the overall effect of the Y-string potential on the spectrum 

does not amount to a mere rescaling of the energy, i.e. of the string tension a in the A-string potential: the Y-string 
potential pushes all of the K=2 states up in energy, some less others more, whereas the A-string shifts some K=2 
states up, others down and leaves third states unchanged relative to the Y-string states. Moreover, the size of the 
energy splitting of the K=2 band states in the Y-string (0.381684) is more than two times bigger than in the rescaled 
A-string potential (0.155752). 

Of course, the ordering of the shells/bands and of states within shells remains identical in both potentials, because 
it is controlled by the permutation symmetry of the states and of the potentials, as well as by the convexity of the 
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[20, 1 + ] 
[70,2 + ] 

[56,2 + ] 
[70,0 + ] 



2A 

~r 

2A 

-4- 



(a) 



[56,0] 



(b) 



(c) 



4.73A 



FIG. 5: Depiction of the energy splitting of the K = 2 states of the hyper-spherical linear potential spectrum due to attractive 
three-body potentials: (a) the first-order perturbation approximation to the first-power expansion of the "three-string" potential 
(see text for definition); (b) exact non-perturbative results of the "complete-string" potential and (c) the (rescaled) A-string 
potential results of J.M. Richard and Taxil 0. The left-hand side of the diagram involving the [20, 1+], [70, 2+], [56, 2 + ], [70, + ] 
multiplets is common to both kinds of potentials, follows the Bowler- Tynemouth rule 2:2:1 to 1%; only the position of the 
[56, + ] multiplet (containing the Roper resonance) is variable. 



TABLE VI: The eigen-energies (in units of y-^JLj ) of the Y-string potential E^. and the A-string potential results of J.M. 
Richard and Taxil [l9l ] rescaled by factor f 2 ^ 3 — 0.843726, see text, for the low-lying K = 0,1,2 states. 



K 


N K 


[SU(6),L P ] 




J C N K ,K,L 








[56, 0+] 


5.1761 


5.1761 


1 





[70,1-] 


6.3160 


6.3160 





1 


[56, 0+] 


7.1360 


7.1360 


2 





[70, + ] 


7.1733 


7.2832 


2 





[56,2+] 


7.2437 


7.3143 


2 





[70, 2+] 


7.3968 


7.3767 


2 





[20, 1+] 


7.5550 


7.4390 



potential in the hyper-radial coordinate. The only qualitative difference between the A- and Y-string potentials that 
we found is that the "Roper" lies lower, as compared with the lowest lying K=2 state, in the A-string than in the 
Y-string potential. 

So we may conclude that, neglecting relativistic and (the very important) HFI effects, the three lowest-lying bands 
of states that form the (only) set of well-established resonances so far, do not allow a clear distinction to be made 
between these two types of potentials. This may, but need not, be a surprise, as these two string potentials have (very) 
different functional forms in configuration space, which we (naively) expected to predict different physics. But it turns 
out that the Bowler- Tynemouth theorem holds for general hyper-radial potentials, not just for the harmonic oscillator 
one, and then the separation of states in the K=2 shell is (tightly) constrained by the permutation symmetry. 

The technical complexities of the Y-string potential, such as the orbital angular momentum dependence, which is 
responsible for some basic features of the spectrum, and the coupling of the radial and angular motions, encountered 
above are difficult to deal with accurately with methods that are not based on the hyper-spherical coordinates. The 
latter technique is not widely familiar to many practitioners in this field, however. The hyper-spherical formalism has 
one unfortunate disadvantage: the relativistic kinetic energy of three quarks is not a function of only the hyper-radius 
and as such the method does not (readily) extend to relativistic energies. For this reason we expect the present results 
to be relevant and applicable only to heavy-quark baryon spectroscopy, such as ccc. 

That may explain the reason why this three-string potential has often been replaced in practical calculations by 
various two-body approximations, that are more easily handled by standard methods. Several two-body approxima- 
tions to this potential have been devised and compared with the exact Y-string potential in various limits [12U13. 22]. 
The results of the two most systematic papers so far [l3], |22j are inconclusive, however, because the former did not 
calculate the spectrum and the latter is based on the hyper-radial approximation, which does not split the various 
SU(6) multiplets within the K=2 shell. 

There is only one possible clue to the shape of the confining potential in the lower end of the baryon resonance 
spectrum, viz. the Roper resonances (abnormally low) mass, that perhaps could be used to draw conclusions about 
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the existence and/or preponderance of one kind of potential over the other. We have shown, however, that the Y- 
shaped string always leads to a Roper resonance that is heavier than the lowest-lying odd-parity resonance, just like 
the A-string. 

This does not mean, however, that there need not be any differences between the spectra of the Y- and the A-string 
potentials, rather, it means that one must go to the higher lying bands of states, and in particular to higher orbital 
angular momentum states, in order to see these differences. It remains to be seen just how high is high enough. 
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APPENDIX A: THE HYPER-SPHERICAL APPROACH 
1. Brief review 

The matrix elements/expectation values of the three-body potential V between the hyper-spherical harmonics 
yK,L, mi {^b) and y K ' tL > t7n '~{P-5) are defined by 



(Vr(R)) lK ][K'] = J y*K,L, m M) V (P> X )yK',L' im ^)dn^ 

= r /2 cos 2 X sin 2 X dx [\m9d9y* KLm ^ 5 )V(p,X)y K , L , m ,(n 5 ). (Al) 
Jo Jo * 



where 



cos 9 — cos dp cos6\ + cos(</?p — ip\) sin9 p sin 9\ 

= " ( i Y ^ P ) ® Y ^x)] L=0 (A2) 

The infinitesimal volume element dV = R 5 dRdfl^ 5 ' and the infinitesimal hyper-spherical solid angle dfl^ are given 
by 

dV = R 5 dRdrt 5 \ (A3) 

dfl {5) = cos 2 X sin 2 X dx sin 9 d9 dfl {3 \ (A4) 

dft (3) = d(/>dn (2 \ (A5) 

dn^ = cos 2 X sm 2 xd X dnW dn { *\ (A6) 

where dO^ 2 ^ is the usual three-dimensional space differential solid angle. Then one has 

J dnW = 47r, J dttW = Sir 2 , J dn^ = n 3 . (A7) 
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The ground state matrix element (Vy)oo(R) is often called the hyper-radial potential Vy(R); the integrals in the 
hyper-radial Y-string potential Vy {R) can be evaluated numerically as 



v Y {R) = {v Y (R)) 00 = J y* fifi (n x )v Y (p,\)y ,oAn P )dn {5) 
= ^3 / Vy(p^)drt 5) , 



8 r r* 

- I cos 2 x sin 2 x dx I Vy (R,x,0) sin 6 d6 

I" JO JO 

<T\/j-R— [ sinOdO [ y/l + \sm9\ sin 2% (cos x sinx) 2 dx 

V 2 TT J Q J Q 

1.2876 <rJ^R. (A8) 



where |sin0| = vl — cos 2 0. One can easily write down the power expansion of the three-body potential Vy; thus 
one obtains an expansion of the potential in (integer) powers of ordinary spherical harmonics and of sin 2%, which 
can be combined into (new) hyper-spherical harmonics yK,L,rm that can be used in Eqs. (|A1[) . (|A8[) by using the 
(hyper-spherical) Clebsch-Gordan coefficients, or simply by brute-force numerical integration. Note that these matrix 
elements are "almost always" diagonal, i.e., that this potential does not connect hyper-spherical harmonics with 
different K, but only of the same K and L, though possibly with different l p and l\ values, due to properties of the 
three-dimensional space spherical harmonics. 

2. K=0,l,2 hyper-spherical harmonics 

Manifestly, in order to be able to evaluate Eq. (|A1[) one must know the explicit form of the hyper-spherical 



harmonics. Simonov, Ref. 17[, has explicitly written down the K = 0,1,2 states' wave functions (hyper-spherical 
harmonics) and examined their permutation symmetryproperties, although the highcr-K symmetrized hyper-spherical 
harmonics have remained a widely unexplored topic [24| . 

The symmetries of the string potential/hamiltonian are: parity, permutation/spatial exchange of quarks, rotation, 
therefore we see that only wave functions with the same P = (— 1)' 1+ ' 2 , L, and symmetry M,S,A can mix with 
each other. Let Py be the ij-th particle permutation/spatial exchange operator. The permutation symmetry can be 
examined using the following properties of the (p, A) set of vectors with mixed symmetry [25[ 

P12P -* -P (A9a) 

F12A — > A (A9b) 

P13P ^i/?-^A (A9c) 

P13A -> \\ (A9d) 

that furnish the basis for the two-dimension irrep of S3. Using the above relations, the second-order scalar, vector, and 
tensor can be constructed and used to obtain the K = 0, 1, 2 hyper-spherical harmonics, see Ref. [171 ]. Thus, it turns 
out that the 5*3 permutation group symmetrized hyper-spherical harmonics correspond precisely to different SU(6)fs 
symmetry multiplets (Young diagrams/tableaux) of the three-quark system. The symmetrized hyper-spherical har- 
monics have been systematically developed only up to K=2, see comments in Ref. [2J]. We show them in Table IViTl 



APPENDIX B: THE HYPER- ANGULAR MATRIX ELEMENTS 
1. First-order expansion of the Y-string potential 



The following is, of course, just a schematic representation of the total matrix elements, as the complete K=2 band 
wave functions may have two components with different partial orbital momentum components (l p ,l\) and identical 
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TABLE VII: The hyper-spherical harmonics, for K = 0,1,2 values. 



K 


n K 


[SU FS (6),L r ] 






hyper-spherical harmonic 





1 


[56,0+] 


(0 x 0) 




yi = ^[Y ®Y ]o 


1 


3 


[70,1-] 
[70,1-] 


(1 x 0)i 
(0 x l)x 




y 2a = 4^/£cosx[Yi® Y ]i 
y-ib = 4y / |sinx[yo(8>Fi]i 


2 


1 


[70,0+] 
[70,0+] 


(0 x 0) 
(1 x 1) 




y 3a = ^(cos X " - anX'Wo ® Yo] 
y-ib = 42= cos y sin v[Yi ® Ydo 


2 


3 


[20, 1+] 


(1 X l)x 




34 = 7= cos x sin x[Vi ® Ki] 1 


2 


5 


[56, 2+] 


(2 x 0) 2 
+(0 x 2) 2 


y 5 = 


-^.-^ [cos 2 x [v 2 s r ] 2 + sin 2 x[v ® y 2 ] 2 ] 


2 


5 


[70, 2+] 
[70, 2+] 


(2 x 0) 2 
-(0 x 2) 2 
(1 x 1) 2 


y 6a = 


[cos 2 X [F 2 ® Yoh - sin 2 X [Yo ® Y 2 ] 2 ] 
y<ib = 7= cos x sin x[Yi ® Yi] 2 



total orbital angular momentum L, such as the K—2 band multiplets [70, + ] and [70, 2 + ] (see above): 

(K(Lf)\ sin2x|sin6»| |i^(Li)) h yp-ang = (i>K(L f )(x)\ sin2x | ^k(U) (x) )h yP (Y Lf \ \ sin (9 \ \ Y u ) ang 

= ( ^K{L f ){x) I sin 2 X I "fix (Li) (x) )h yP 

x ((*!/ ® Mf/ I I sin^l I (Zli ® « 2i )^)ang. (Bl) 

We may drop the indices for the initial and the final state total orbital angular momenta L j = Li = L; Mf = M{ = M, 
as the ang ular matrix elements are diagonal in those quantum numbers. We use the spatial wave functions available in 
Refs. [l6l IT7I l2lT| to calculate the hyper-angular matrix elements (i/'k | )hyp-ang = (V'K I sin2x | ipn ) hyp- ang 
and take the angular matrix elements (Yj, | | sin^l | Yl )ang from Table Hfto calculate the total hyper-angular matrix 
element (ipK \ | sin 2% | sin#| | ipx )hyp-ang- We list the radial-, angular- and total matrix elements' values, as well as 
the results of some intermediate steps, in Table [VTlTl 



and 



TABLE VIII: The values of the three-body potential hyper-angular matrix elements (-Sr)hyp = ( sm 2x)h yP 
(sin 2x I sin ^|)hy P -ang, for the K = 0,1,2 states (in all the partial S-, P- and D-waves). In the [56, 2 + ] and the [70, 2 + ] en- 
tries "diag." denotes diagonal terms of the type (2 x 0) 2 x (2 x 0) 2 and (0 x 2) 2 x (0 x 2) 2 , whereas "off-diag." denotes 
off-diagonal terms (2 x 0) 2 x (0 x 2) 2 and (0 x 2) 2 x (2 x 0) 2 . These two kinds of contribution must be separated because they 
have different angular matrix elements ^(| sin#|). 



K [SU FS (6),L^] (l p x I x )l N 



[56, H 
[70, V 
[56, H 
[20, 1 H 



[70, H 
[70, 0" 



[5^ 

[56, 2+ 
[56, 2+ 



[70, 
[70, 

[70, 

[70, 



(0 x 0)o 
(1 x 0)l 
(0 x 0)o 
(lxl)i 



^(sin2x) ^(\sm6\) (sin 2 X | sin gj) = 



(0 x 0)o 
(1 x l)o 



diag. 
off-diag. 

(2 x 0)2 
+(0 x 2)2 



75tt 



ill 

1 

2 r, 



diag. 
off-diag. 

(2 x 0)2 
-(0 x 2)2 
(1 x 1)2 



gg 

75tt 



1G 



iiii 

41) 



75 
44 

75 



1 
75 

r>2 
75 
52 



0.6667 
0.6667 
0.6667 
0.8000 



0.5333 
0.5333 



0.6400 
-0.0533 

0.5867 



0.6400 
0.0533 

0.6933 

0.6933 
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2. The Y-string potential 



The total matrix elements may have two components with different partial orbital momenta (l p ,l\) and identical 
total orbital angular momentum L, e.g. the K—2 band multiplets [70, + ] and [70,2+]: where 

(^(L / )(x)lv / l + sin2 X |sin0||^ (ii) ( X )) hyp = Qsin2 X ) \<Pk,l(x)\ Vl + sin2 X | sin0| d X - (B2) 

and we are left with a double integral over the product of even-L order Legendre polynomials Pl{x), the potential 
•v/l-j- sin2 X | sin(9| and the hyper-angular wave function squared |V'-R",l(x)| 2 : 

iy dcos9P L (cos9) (^K(L f ) I V 1 + sin2 X |sin6>| \ ^ K {L i ) ) = \J dcos9P L (cos9) J d X Qsin2 X ^ 

x |V^,L(x)lVl + sin2 X |sin0|. (B3) 

We use the K=0,l,2 hyper-angular wave functions available in Appendix I A 21 to calculate these hyper-angular matrix 
elements and take the angular momentum coefficients from Table HI The integral can be rewritten in terms of new 
variables z — cos 2 X and x = cos 9 as follows 

~J d cos 9P L (cos 9) (ip K (Lf)(x) I V 1 + sin2 X | sin0| \ip K (L i )(x))h yP = \J dxP L (x) J \J\- z 2 dz 

x \1> K ,L Wl 2 \/l + Vl- z 2 y/T^, (B4) 

and evaluated numerically; we list the results in Table HXl 



TABLE IX: The values of the three-body potential hyper-angular matrix elements (yl + sin 2\\ sin #|)hy P — ang. for the K = 
0,1,2 states (in all the partial S-, P- and D-waves). 



K [SU FS (6),L P ] 


(a/1 + sin 2 X sin<9 ) hyp _ ang . 


[56,0+] 


1.2876 


1 [70,1-] 


1.2876 


[56, 0+] 


1.2876 


2 [20, 1+] 


1.3402 


2 [70,0+] 


1.2350 


2 [56, 2+] 


1.2560 


2 [70, 2+] 


1.2981 



3. The complete string potential 

The boundary Eqs. ([50]) can be recast in the new variables z — cos2 X and x — cos 9 of Fabre de la Ripelle and 
Lassaut, Ref. [22j, that are particularly useful in the integration over the solid hyper-angle, see Fig. [6) 

, . 5x 2 -8x 4 + V3xVl-x 2 

zi(x) = 



z 2 (x) 



1 - 3 x 2 + 8 x 4 + V3 x Vl ~ x 2 ' 

-5x 2 + 8x 4 + VZxVl -x 2 
-1 + 3a; 2 - 8.x 4 + V3 x v 7 ! - x 2 ' 



, , 2 + x 2 + V4- 5x 2 + x 4 

Z 3 (x) = ; (B5) 

V ' -7 + x 2 + V4-5x 2 + x 4 
Due to the symmetry of the integrand we may take twice the integral over x from to 1 ; moreover the integral can 
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Z=C0S2;tf 



1 

0.5 

0.5 

-1 





X=COSr? 



-1 



0.5 



0.5 



FIG. 6: The boundary in the z = cos 2\ vs. x — cos 9 plane, between the regions in which the two- and the three-string 
potentials are appropriate, see Eqs. (|B5|) . 



be rewritten in terms of the new variables z = cos 2\ and x — cos 9 as follows 



d cos 6P L (cos 6) (ipK(L f )(x) I V(x,0) | ip K {Li)(x) >h yP 



pi pi pzi(x) y 

dx I dz + dx dz 

n Jz 3 (x) Ji Jz 3 (x) j 



'l-x 2 



x (^|7p^ + yi7|^T=^ 

+ Q) ^ P L {x)dx \i> K{L) {z)\ 2 ^/l-zHz 

x (71^+^1 + ^-^3(1-^)), (B6) 
where 21,3(2;) are given by Eqs. (|B5[) defining the boundary in the right-hand-side half of the x vs. z plane. 

TABLE X: The values of the three-body potential hyper-angular matrix elements (Ktring)hyp-ang. for the K = 0,1,2 states (in 
all the partial S-, P- and D- waves). 



K [SU FS (G),L F ] (Ktring)hyp- 






[56,0+] 


1.289126 


1 


[70, 1"] 


1.289126 





[56, 0+] 


1.289126 


2 


[20, 1+] 


1.340371 


2 


[70, 0+] 


1.240090 


2 


[56,2+] 


1.258379 


2 


[70, 2+] 


1.298492 
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K 


a 


Oi 


a2 





3.82 


5.26 


6.54 


1 


4.66 


5.99 


7.19 


2 


5.43 


6.67 


7.81 



TABLE XI: Energy eigenvalues of the three-quark states in the linear hyper-radial potential for the ground- (do), the first- (ai) 
and the second (0,2) radially excited states, in natural units (see text). 



APPENDIX C: THE HYPER-RADIAL SCHRODINGER EQUATION 
1. Solving the hyper-radial equation 

The eigenvalues of the three quark states can be obtained by solving the eigenvalue problem of the following equation 



d 2 5 d K(K + 4) [3 

d^ + RdR-^B^ +2 ^ E - 2 ^r R 



1p[K] = 



By changing scale of the hyper-radial coordinate one can write 

d 2 5 d K(K + 4) 



dx 2 x dx 



+ a — x 



l/>[K] = 



with 



R = ax 



-1/3 



a = 2(iEa 2 



(CI) 

(C2) 

(C3) 
(C4) 
(C5) 



Regular solution can be written ip[K] = x u [k]( x ) with U[k](0) — const. The equation for u<k] can be written as 

u [K ] = (C6) 



d 2 2K + 5 d 

■7-0 ^ r~ + a - x 

dx z x dx 



(C7) 
(C8) 



The above equation is solved with the boundary condition 

u(Q) — finite 
ti(oo) — * 

The obtained eigenvalues are tabulated in Table IXII 

2. Evaluation of the potential's hyper-radial matrix element 

We do not need the hyper-radial wave functions in order to complete the calculation of the radial matrix elements; 
rather, we use the virial theorem 



V(k) 



k + 2 



E{k) 



to determine the expectation value (V(k = 1)) = §-E(fc = 1), where k is the power of r in the potential V(k) ~ r k , in 
terms of the energy eigenvalue E{k = 1) of the unperturbed hamiltonian. 
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"the Y-string and the A-string potentials are essentially indistinguishable" 

It would be ideal if these two kinds of string potentials were to predict (dramatically) different baryon spectra that can be 
easily differentiated by experiment, but that turns out not to be the case 

It has been claimed that this approximation is exact [%|, however, in the case(s) of zero individual orbital angular momenta 
(l p , i.e. for the radial excitations of the S-wave ground state, and presumably a good one for low values of the orbital 
angular momenta. We shall show below that this claim is incorrect, however, 
but not the N=0,1 shells states. 



